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ABSTRACT 

We present the first mid-infrared interferometric measurements of FU Orio- 
nis. We clearly resolve structures that are best explained with an optically thick 
accretion disk. A simple accretion disk model fits the observed SED and visibili- 
ties reasonably well and does not require the presence of any additional structure 
such as a dusty envelope. The inclination and also the position angle of the disk 
can be constrained from the multibaseline interferometric observations. Our disk 
model is in general agreement with most published near-infrared interferometric 
measurements. From the shape and strength of the 8-13/xm spectrum the dust 
composition of the accretion disk is derived for the first time. We conclude that 
most dust particles are amorphous and already much larger than those typically 
observed in the ISM. Although the high accretion rate of the system provides 
both, high temperatures out to large radii and an effective transport mechanism 
to distribute crystalline grains, we do not see any evidence for crystalline silicates 
neither in the total spectrum nor in the correlated flux spectra from the inner 
disk regions. Possible reasons for this non-detection are mentioned. All results 
are discussed in context with other high-spatial resolution observations of FU Ori 
and other FU Ori objects. We also address the question whether FU Ori is in a 
younger evolutionary stage than a classical TTauri star. 

Subject headings: Accretion Disks - Techniques: Interferometric - Circumstellar 
Matter - Stars: Formation - Stars: Pre-Main Sequence - Stars (Individual): FU 
Orionis 



1 Based on observations made with ESO Telescopes at the Paranal Observatories under program ID 
074.C-0209 and 274.C-5032 
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1. Introduction 

FU Orionis is the prototype of a small, but quite remarkable class of low-mass Young 
Stellar Objects (YSOs) normally referred to as FU Ori objects (FUORs). For the first mem- 
bers of this class an outburst in optical light of up to 4-6 magnitudes over short time scales, 
followed by a decrease in luminosity over several years or decades, was observed (Herbig 
1977). Other objects were included in the class as they shared common specific spectroscopic 
features, e.g., double-peaked line profiles, a spectral type varying with wavelength and often 
CO bandhead absorption features. As for one object (V1057 Cyg) the pre-outburst spectrum 
is known and as it resembles that of a classical TTauri star (Welin 1971), it is commonly as- 
sumed that FUORs should be low-mass YSOs. Most observational data can be explained by 
the presence of an accretion disk surrounding the young stars (Herbig et al. 2003, however, 
found that some spectral properties are well explained in the context of rapidly rotating 
late- type stars with strong stellar winds). A dramatic temporal increase in the accretion 
rate, where the disk outshines the star by several orders of magnitude, can account for the 
observed outbursts in luminosity. Several scenarios, possibly triggering such an increased 
accretion rate, have thus far been proposed. They include (a) interactions of binary or mul- 
tiple systems where tidal forces disturb the circumstellar disk (Bonnell & Bastien 1992), (b) 
planet-disk interactions, where thermal instabilities in the disk are caused by the presence 
of a massive planet (Lodato & Clarke 2004), or (c) thermal instabilities in the disk alone 
(Bell et al. 1995). For a detailed overview concerning the FU Ori phenomenon we refer to 
Hartmann & Kenyon (1996). Apart from revealing the mechanism leading to the observed 
outbursts, it is also important to investigate whether all TTauri stars undergo such epochs 
of enhanced accretion or whether FUORs are a special class of YSOs. Most observations 
of classical TTauri stars show that the derived accretion rates of 10~ 10 — 10~ 7 M Q /yr (e.g., 
Gullbring et al. 1998) are not sufficient to build up a low-mass star over time scales of a few 
Myr. Even if part of the matter is supposed to be accreted in the very early phases of a 
YSO, FUOR-phases might provide an elegant solution to this problem as they would speed 
up the accretion process. 

As so far only a dozen or so well established FUORs are known and as the FUOR- 
phase has significant impact on the young star-disk system it is crucial to combine as much 
observational information as possible for these objects to eventually derive a coherent theo- 
retical picture. Based on near-infrared (NIR) and/or mid-infrared (MIR) interferometry as 
a technique to study the inner few AU of the accretion disks new insights to some of the best 
studied FUORs were provided recently. Millan-Gabet et al. (2006) found that accretion disks 
alone can not reproduce the SED and observed low K-band visibilities for V1057 Cyg, V1515 
Cyg and Z CMa-SE simultaneously. They concluded that additional uncorrelated flux may 
arise due to scattering by large dusty envelopes. Abraham et al. (2006) presented the first 
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VLTI/ MIDI observations of VI 647 Ori whose eruptive behavior suggests that it is either an 
FUOR or an EX Lupi (EXor) type object. In this case it was possible to fit both, the SED 
and the observed MIR visibility with a simple disk model with moderate disk flaring. For 
FU Ori itself it was Malbet et al. (2005) who analyzed a wealth of NIR interferometric data. 
They showed that the NIR visibilities and the SED could be fitted with two models: One 
consisting simply of an optically thick and geometrically thin accretion disk and a second 
one consisting of an accretion disk and an embedded "hot spot". From their error statistics 
these authors concluded that the latter model was more likely. 

In summary it becomes clear that up to now no coherent picture can be derived from 
the interferometric observations of FUORs and the group of objects seems to be rather 
inhomogeneous. 

In this paper we present the first MIR interferometric measurements of FU Orionis. The 
data are thus complementary to the NIR observations of Malbet et al. (2005). In section 2 
we briefly describe the observations and the data reduction process. In section 3 we discuss 
the findings derived from the N-band acquisition images. The 8-13/zm spectrum of FU Ori 
is analyzed in section 4. In sections 5 and 6 we discuss the results from the interferometric 
measurements, i.e. the visibilities and the correlated flux spectra, respectively. A simple 
analytical disk model and its fit to the SED and the visibilities is presented in section 7. 
Finally, we summarize our conclusions and mention some future prospects in section 8. 

2. Observations and Data Reduction 

The observations were carried out between October 31 st and November 4 th 2004 with 
the Mid-Infrared Interferometric Instrument (MIDI) at ESO's Very Large Telescope Inter- 
ferometer (VLTI) on Paranal/Chile. Together with the Keck Interferometer Nuller that 
recently started to produce first scientific results (Mennesson et al. 2005), MIDI is currently 
worldwide the only instrument able to conduct spectrally resolved interferometric obser- 
vations in the mid-infrared (Leinert et al. 2004). For the observations MIDI was used in 
high-sens mode using the NaCl prism as dispersive element yielding a spectral resolution 
of 30. The maximum projected baseline was 86.25m and the minimum projected baseline 
44.56m leading to an angular resolution at 10/xm of 0."029 and 0."056, respectively. For an 
assumed distance of 450 pc these values correspond to 13.1 AU and 25.2 AU. A journal of 
the observations including projected baselines, position angles and calibrator stars is given in 
Table 1. For completeness we also mention the observations from December 2004 although 
they are disregarded in the following sections. These observations had almost exactly the 
same baseline and position angle as the October observations but showed in general a lower 
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level of total and correlated flux due to rather poor observing conditions in terms of seeing 
and transparency. The calibrator stars for the data reduction were chosen by analyzing all 
calibrator stars observed over the whole night. We selected those showing a good agreement 
in their transfer functions, i.e. their instrumental visibilities after the assumed sizes were 
taken into account. 

The data reduction was carried out with the software package MIA+EWS-1.3. This 
software consists of two independent reduction programs (MIA and EWS) which were both 
applied for comparison. The software, further information and manuals can be downloaded 
from the Internet 2 . A general description of the basic data reduction steps is also given in 
Leinert et al. (2004) and we refer the reader to this paper. As the results derived with MIA 
and EWS agreed quite well (overall differences in the correlated flux < 8%), we decided to 
show only plots resulting from the MIA package. 



3. The MIDI Acquisition Image 

FU Ori was found to be a binary system by Wang et al. (2004). Reipurth & Aspin 
(2004) confirmed the detected companion and concluded that it was a young star of spectral 
type K showing clear NIR excess. In the same paper it was also stated that it was very 
unlikely that the binary component (FU Ori S) triggered the outburst of FU Ori observed 
in the late 1930s. 

Knowing about the existence of the fainter companion the integration time of some MIDI 
acquisition images was increased in order to derive N-band photometry for both components. 
In three images FU Ori S was clearly visible (Figure 1) and aperture photometry could be 
applied to the observations. The results are summarized in Table 2. Interestingly, FU Ori S 
shows relatively a higher N-band excess than FU Ori itself. This can be explained in terms 
of differences in the geometry of an assumed circumstellar disks (e.g., larger flaring angle or 
smaller disk inclination). We furthermore derive a separation between the components of 
0."484±0."01 and a position angle of 162.5°±4.1° (measured from north eastwards) which is 
in good agreement with the results from Wang et al. (2004) and Reipurth & Aspin (2004). 
Our errors are standard deviations based on Gaussian fitting of 3 independent images. For 
an assumed distance of 450 pc the separation corresponds to 217.8±4.5 AU. 



2 http://www.strw.leidcnuniv.nl/^nevcc/MIDI/index.html 
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4. The Total Uncorrelated MIR Spectrum 

4.1. Comparison to Spitzer observations 

Figure 2 shows the MIR spectrum of FU Ori. The spectrum obtained with MIDI is 
the average over the first three observing nights using both Unit Telescopes (UTs). For 
comparison Spitzer IRS data are plotted additionally. Those data were available from the 
Spitzer Data Archive. For the data reduction we used the droopres intermediate data product 
processed through the SSC pipeline S12.0.2. The SMART reduction package developed by 
the IRS Instrument Team at Cornell (Higdon et al. 2004) was used to extract the spectrum. 
Within the error bars the spectra agree quite well. A broad and weak silicate emission 
feature is present between 8-13/im However, the averaged MIDI spectrum appears to be a 
little flatter than the more roundish Spitzer spectrum. 



4.2. Possible dust composition 

In contrast to other YSOs the observed silicate emission feature of FU Ori is rather weak 
(see also Hanner et al. 1998). Such a broad and flat MIR dust spectrum can be explained by 
grains that already underwent some coagulation process (Bouwman et al. 2001). To test this 
hypothesis we derived a possible dust composition by fitting a dust model to the normalized 
and continuum subtracted Spitzer spectrum from 6-13/xm. This analysis method is similar 
to approaches that were successfully applied in previous studies for determining the dust 
composition of observed dust emission features of other YSOs (e.g., Bouwman et al. 2001). 
In the model we assume that the observed emission can be computed from the sum of the 
emission of individual dust species. The species we used are summarized in Table 3 together 
with references for their optical properties. We furthermore took into account different grain 
sizes as for each of the dust species the opacities were calculated for 0.1, 1.5, and 6.0/im 
sized particles. To estimate the possible contribution of Polycyclic Aromatic Hydrocarbons 
(PAHs) to the emission we used a template spectrum based on profiles that were derived 
from observations (Peeters et al. 2002; van Diedenhoven 2004). Finally, in order to reproduce 
the observed spectrum we fitted the following emission model to the Spitzer data using a 
least square minimization 

F v - B u {T cont )C, = B„(T Awt ) (YsYI CijitA + C PA h4 PAH • (1) 

V i=i j=i / 

Here, F u is the observed flux, B„(T cont ) denotes the Planck function at the temperature of 
the underlying continuum, B v (Td U st) the Planck function at the temperature of the silicate 
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grains, k 1 ^ is the mass absorption coefficient for the silicate species j with the grain size 
i, I^ AH is the template emission spectrum for the PAHs , and C , Cij, and Cpah are 
the weighting factors for the continuum, the silicate emission and the PAH contribution, 
respectively. A more detailed description of the dust model and plots showing the different 
dust opacities is given in Bouwman et al. (2006). 

The model results are shown in Figure 3 and Table 4 summarizes the corresponding 
dust composition. As for the underlying continuum we could fit a temperature of 880 K, we 
believe that this optically thick emission most likely arises from the inner parts of the disk 
close to the star. The outer and thus cooler regions of the disk would then be responsible for 
the weak silicate emission for which we fitted a temperature of 230 K to the optically thin 
emission layer. The results strongly support the idea that most particles are amorphous and 
much larger in size than those in the ISM. PAHs, silica particles and crystalline silicates do 
not seem to be present at all. Although the derived fractions of different dust species do 
depend on the model applied to the data (e.g., in terms of the precise grain sizes considered 
and also in terms of particle structures) the main conclusions, that we do not find evidence 
for crystalline silicates and that the grains are considerably larger than those found in the 
ISM, remain unaltered. For a more detailed discussion about the influence of different model 
parameters we refer to Voshchinnikov et al. (2006) and van Boekel et al. (2004). 

As the disk is much brighter in the optical and NIR than the central star (section 7), 
stellar radiation is probably not able to produce the observed weak emission feature in an 
optically thin disk surface layer. However, it was shown that the silicate feature of FU Ori 
can be reproduced taking into account self-irradiation of the disk (Lachaume 2004). In this 
case the hot inner parts of the accretion disk serve as flux source and illuminate the disk 
surface layer under a certain angle. 

5. The Visibilities 

In Figure 4 and Table 5 the visibilities measured for FU Ori at three different baselines 
are shown. The errors are computed as standard deviations resulting from the use of different 
calibrators for each night (see Table 1). In most of the following plots the position of the 
atmospheric ozone band is indicated as this part of the spectrum suffered sometimes from 
imperfect corrections during the data reduction process. Here, small "dips" and "bumps" in 
otherwise flat spectra are remnants from the data reduction. 
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5.1. Qualitative assessment 

The calibrated visibility of the shortest baseline (UT2-UT3) remains almost constant 
at ^0.95 over the whole wavelength range and the source is only marginally resolved. As 
expected, for the longest baseline (UT2-UT4) the lowest visibility is observed ranging from 
«i0.73 at 8.3;um to ^0.6 at 13/im. For the intermediate baseline (UT3-UT4) the visibility 
shows a slight increase from f»0.73 at 8.3yum to ?^0.85 at I3fim. The fact that the object 
is clearly resolved with two baselines in the MIR supports the argument of Malbet et al. 
(2005) that we are observing an extended circumstellar structure, i.e. a disk, and not only 
the stellar photosphere (Herbig et al. 2003). 

For the two baselines UT2-UT3 and UT2-UT4 the observations are consistent with ex- 
pectations from thermal disk emission as despite the decreasing resolution for longer wave- 
lengths the visibilities indicate larger sizes for the emitting regions (see also section 5.3.). 
The observed increase in visibility for the UT3-UT4 baseline implies that for this baseline, 
however, the object appears smaller at 13/xm than at 8/im. This seems difficult to imagine in 
the context of a circumstellar disk as emitting source and certainly a second observation for 
this baseline configuration seems eligible. Theoretically it is possible that the photometric 
measurements carried out directly after the interferometric observations are corrupted due to 
technical problems or different weather conditions. This in turn might then lead to a change 
in the visibility function. However, the calibrators observed at this night did not show any 
sign of poor photometric measurements and their transfer function was very stable over the 
whole night. In addition, also other YSOs with circumstellar disks showed an increasing 
visibility for certain baselines when observed with MIDI. Thus, since we do not find any 
evidence for excluding this dataset due to bad quality we decided to keep it in our analyses. 



5.2. Comparison to MIR visibilities of HAeBe stars 

As thus far no MIDI visibilities for TTauri stars have been published we are limited to a 
comparison between FU Ori and circumstellar disks around Herbig Ae/Be stars (HAeBes). 
It shows that these stars as well as models applied to them normally show a prominent 
drop in the visibility between 8/xm and 10/zm from where the curve remains almost constant 
(Leinert et al. 2004). Qualitatively, this drop results from the intensity distribution of the 
passive disks around the HAeBes: At the short wavelength end the hot inner rim of the 
disks provides an overproportional contribution to the flux and is at the same time confined 
to a small spatial region leading to a high visibility. Most of the rest of the MIR emission 
originates from a large area of the hot, illuminated surface layer of the flared circumstellar 
disk. The visibilities for FU Ori show hardly any wavelength dependence and are very flat 
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from 8-13/xm regardless of baseline. Thus, its seems as if the flux distribution of FU Ori 
is smoother and the visibilities show no sign of a significant contribution from a hot inner 
rim. However, as FU Ori is surrounded by a heavily active accretion disk where the disk 
alone produces the majority of the observed flux at almost all wavelengths (see section 7) 
differences in the intensity distribution and thus in the visibility can be expected. 

Apart from differences in the shape of the visibility curves most HAeBes are much better 
resolved, i.e. show lower visibilities than FU Ori (Leinert et al. 2004). This, however, can at 
least partly be explained by the distance to these objects which is in general less than 200 
pc. FU Ori on the other hand has an assumed distance of 450 pc, and if it was closer to the 
Earth we would observe lower visibilities also for this object. 



To derive a simple model for the geometry of the emitting regions we assume a simple 
Gaussian brightness distribution for each baseline. This is a reasonable first approxima- 
tion for objects showing high visibilities. The FWHM, and hence the physical size, of this 
Gaussian in arcsec can be computed by 



derived from the projected baseline B and the wavelength A. Equation 2 results from a 
simple Fourier transformation of the assumed brightness distribution. We computed the 
FWHM for all three baselines at three different wavelengths (9.0, 11.0, and 12.5/xm). For 
this we averaged for each baseline 5 visibility points from Table 5 centered on the specified 
wavelengths. Table 6 gives the resulting sizes of the emitting regions in AU for an assumed 
distance to FU Ori of 450 pc. The results are also visualized in Figure 5 where the FWHM 
are shown in their orientation on the sky. As expected for thermal disk emission the FWHM 
increases with wavelength for a given baseline. The only exception is the 12.5/im size for the 
UT3-UT4 baseline which is surprisingly a little smaller than that seen at 11.0/iin. 

In addition to size estimations for each individual baseline our measurements based on 
different position angles allow to constrain the geometry of the disk. We fitted the derived 
FWHM with an ellipse for each of the considered wavelengths in order to derive a simple 
model for the spatial orientation of an assumed disk-like structure. The resulting best fit 
ellipses are overplotted in Figure 5 and their parameters are summarized in the lower half of 
Table 6. As one would expect from a disk-like structure showing a decrease in temperature 



5.3. Geometry of the emitting regions 




(2) 



where V(f) is the measured visibility for a certain spatial frequency f — ^ (in arcsec : ) 
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with radius the semimajor axes of the ellipses increase with wavelength. However, due to 
the increasing visibility of the UT3-UT4 baseline the semiminor axis and also the total area 
of the ellipse are slightly smaller at 12.5/xm than at 11.0/zm. 

An inclination angle for an assumed circular disk can be derived by computing the 
arccos of the ratio of the semiminor and semimajor axis. The resulting angles agree very 
well (55.4°±2.4°) and are also in good agreement to what was found for FU Ori disk models 
based on NIR observations (Malbet et al. 2005). By combining the derived inclination 
angle with rotational velocity measurements in the optical spectral region it is possible to 
estimate the central mass. Based on Kenyon, Hartmann & Hewett (1988) we derive for FU 
Ori M/Mq r^O.36 which is reasonably consistent with typical values for low-mass pre-main 
sequence stars. 

In addition to the inclination the fitted ellipses provide also information on the position 
angle of the disk. The position angles (measured from north eastwards) we find for the 
semimajor axes at different wavelengths show a larger scatter (109. Fill. 6°) than the incli- 
nations and differ significantly from the NIR findings of Malbet et al. (2005) (see also section 
7). It seems as if the position angle becomes smaller for longer wavelengths, although one has 
to keep in mind that, again, at least for the UT3-UT4 baseline the increasing visibility is the 
main reason for the observed rotation. At this point we leave it to future investigations to 
derive a 3-D disk model based on a more complex dust distribution possibly able to confirm 
the apparent changes in the position angle for different wavelength regimes. 



6. The Correlated MIR Spectra 

6.1. The origin of the correlated flux 

The correlated flux is directly linked to the total flux (i.e. the flux from a single UT 
telescope) via the visibility: 

F corr (\) = V(\)-F total (\) (3) 

Figure 6 depicts the results for the correlated flux for the three different baselines. The total 
spectrum is plotted for comparison. 

Taking into account the spatial resolution of the three baselines it is possible to estimate 
from where the correlated 8-13/xm flux originates. Our UT2-UT3 baseline has a spatial 
resolution of ~25 AU at the distance of FU Ori, and if we assume an average visibility of 
0.95 we hence know that 95% of the 8//m-13//m flux must come from within this 25 AU. 
Similarly, we find that for the UT2-UT4 baseline 65% of the flux is emitted from the inner 
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13 AU. For the UT3-UT4 baseline 75% of the 8-13^m flux comes from the inner 20 AU. 



6.2. Lacking crystalline silicates? 

The spatially resolved MIR spectra of the MIDI instrument make it possible to study 
the radial dependence of the dust composition in the protoplanetary disk (van Boekel et 
al. 2004). In section 4 we showed that the spatially unresolved (total) spectrum of FU Ori 
can be fitted with mainly large amorphous dust grains. The shape of the spatially resolved 
(correlated) flux spectra contains information about the dust composition of the inner parts 
of the disk. 

To compare the shape of the correlated spectra to the total spectrum we first subtracted 
the continuum for which we fitted a straight line between the flux at 8.25 and 12.95/im. 
Then we normalized all continuum subtracted correlated flux spectra and the total Spitzer 
spectrum by 

F = Fsub (A) 

± norm „ V / 

f mean, sub 

where F su ^ and F mean ^ u i, denote the continuum subtracted spectra and their mean value, 
respectively. By applying this formula all correlated spectra conserve their shape and can 
be compared to the total spectrum (Figure 7). To identify those parts where the spectra 
differ significantly we furthermore computed the deviation from the total Spitzer spectrum 
for each wavelength in units of the data error: 

Fnorm,corr Fnorm,total / c \ 

Vdev = (OJ 

<7 corr 

F nor m,corr and F norm ^otai are the normalized correlated spectra and the normalized total 
spectrum and u corT denotes the standard deviation of the correlated spectra (MIDI errors) 
at a given wavelength. We applied the same normalization factors to the errors as we did to 
the spectra. Figure 7 shows that no clear deviations from the total unresolved spectrum are 
observed for any of the correlated spectra. This means that with the given spatial resolution 
and sensitivity no significant chemical difference in the dust composition is observed and the 
spectra of the inner parts of the disk look very similar to the total unresolved disk spectrum. 

This finding is somewhat surprising. As seen in section 4 the silicate dust grains appar- 
ently already underwent noticeable coagulation. Thus, the disk does not consist of purely 
ISM dust anymore and must already have a significant age. In addition it is known that 
apart from grain growth also thermal annealing takes place within protoplanetary disks and 
transform the amorphous silicates into crystalline silicates (see, Henning et al. 2006; van 



- 11 - 



den Ancker 2000). However, neither the total spectrum nor the correlated spectra indicate 
the presence of crystalline silicates. This is even more surprising considering the high tem- 
peratures in the disk interior due to viscous heating of this actively accreting star (section 
7). Annealing processes transforming ISM amorphous dust into crystalline dust set in in 
disk regions where T > 800K (e.g., Gail 2004). In addition to that, high accretion rates 
directly relate to a faster vertical and horizontal mixing within the protoplanetary disk (Gail 
2001). This effect also eases the detection of crystalline silicates within the accretion disk 
surrounding FU Ori. But why don't we find any evidence for their existence? 

A possible explanation might be that they simply cannot be detected as the related 
feature is rather weak and simply dominated by the underlying continuum. This hypothesis 
can be checked by analyzing the longer wavelength regime of the Spitzer spectrum as shown 
in Figure 8. As at longer wavelengths the temperature of the continuum is decreasing 
the contrast between the continuum and possible emission features due to the presence of 
crystalline silicate particles increases. Hence, spectral signatures of these particles should 
be more easily detectable. Figure 8, however, shows that also between 14 and 34/im the 
spectrum of FU Ori does not show any clear signs of crystalline silicates. Normally crystalline 
forsterite particles have emission bands at 16.4, 23.9, 27.7 and 33.8/im and enstatite grains 
show features at 18.8, 21.5, and 24.5/iin (Molster & Kemper 2005). If one argued that 
possibly around 27.7-27.8/zm a weak forsterite feature was present one should keep in mind 
that the features at 23.9 and 33.8/im are normally stronger (Min et al. 2005) but are not 
seen at all. Thus, also the longer wavelength MIR spectrum of FU Ori shows that crystalline 
silicates are not present. This shows that apparently the non-detection is at least not only 
due to contrast problems between the emission features and the underlying continuum. 

What else can prevent us from detecting crystalline silicates? In case FU Ori underwent 
its first (of possible multiple) outburst (s) and phase (s) of enhanced accretion, it is possible 
that not enough time has yet elapsed for the crystalline silicate particles to reach the disk 
surface from where they can be detected by means of MIR spectroscopy. This hypothesis is, 
however, questionable as the large sizes of the amorphous silicates indicate that enough time, 
at least for dust coagulation, has elapsed. Recently, Dullemond, Apai & Walch (2006) found 
that assuming the crystallization by thermal annealing happens at the very early phases of 
disk evolution, the level of crystallinity is linked to the properties of the molecular cloud core 
from which the disk formed. It was found that rapidly rotating cores produce rather massive 
disks with a high accretion rate but a low level of crystallinity. As we observe these properties 
for FU Ori the proposed model, although yet very simple and with limited predictive power, 
might provide an interesting way to explain at least partly the lack of crystalline silicates in 
the FU Ori accretion disk. Finally, another possibility might be that the vertical and radial 
mixing within the disk does not work as efficient as expected. In consequence not enough 
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crystalline dust grains reach the disk surface at radii where also the contrast between the 
emission feature and the dust continuum allows their detection. 



6.3. Comparison to TTauri stars and HAeBe Stars 

The results from the correlated flux allow a comparison to disks around other YSOs that 
have also been observed with VLTI/MIDI. Not only seem crystalline silicates be commonly 
present in protoplanetary disks, but for three intermediate-mass HAeBes stars van Boekel 
et al. (2004) found that the degree of crystallization within the disk depends on the distance 
from the central star. The inner parts of the disks are more crystalline than the outer 
disk regions. This can naturally be explained by assuming that in the inner disk regions 
where the disk has higher temperatures dust processing and annealing occurs and that over 
time crystalline particles are transported outwards. Not only for intermediate mass stars 
but also for the TTauri star RY Tau the same radial dependence of the crystallinity was 
observed (Schegerer et al. 2006). Thus, crystalline silicates should in general be more easily 
detectable in the correlated spectra of protoplanetary disks. For reasons mentioned in the 
previous section it seems thus surprising that the disk around FU Orionis does not seem to 
contain any detectable amount of crystalline dust. 

As described in the introduction FUORs are believed to be low-mass or TTauri-like stars 
in a state of enhanced accretion. From sub-millimeter and millimeter observations Weintraub 
et al. (1991) concluded that FUORs are possibly in an evolutionary youger age than classical 
TTauri stars. This younger age is also supported by theoretical models where the protostellar 
disk is continuously replenished by a dusty envelope. This material infall on the disk can lead 
to thermal instabilities in the disk and thus provides a possible mechanism for explaining 
the observed outbursts Bell et al. (1995). The puzzling thing is now that apparently in the 
disk surrounding FU Ori the amorphous grains have already grown considerably. Hence, 
the disk can not be in its native dust composition and might already have a significant 
age. Furthermore, from the Spitzer FEPS (Formation and Evolution of Planetary Systems) 
legacy program it is known that apparently the amount of crystalline silicates in TTauri disks 
is not correlated with the average size of amorphous dust particles (Bouwman et al. 2006). 
Hence, crystalline dust particles are created in the early phases of disk evolution where the 
accretion rate and the disk temperature are thought to be higher. If thus the FUOR-phase 
is indeed common to most/all TTauri stars then the non-detection of crystalline particles 
cannot be explained in the current picture of protoplanetary disk evolution. 
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7. A Simple Disk Model 



7.1. 



The SED 



The present paper is more intended to show the observational results rather than to 
present a sophisticated disk model possibly able to match the observed SED and visibilities 
with highest accuracy. Thus, instead of applying a complex radiative transfer model we 
decided to use in a first step a simple analytical disk model. In this model the effective 
temperature and also the surface density are prescribed by the following broken power law 
distributions 



For the temperature law the temperature at 1 AU and the exponent q are free parame- 
ters. For the surface density the total disk mass M D and the exponent p are free parameters. 
As the integrated surface density must equal the total mass of the disk the radial dependence 
of the surface density E(r) can be computed. The shape of the total SED (and finally also 
the observed visibilities) require a double power-law to fit the data. At a radius of 3 AU the 
temperature and surface density distribution change. However, we made sure that for both 
parameters the transition was smooth and continuous. Figure 9 shows the SED and the 
resulting fit from the simple disk model. The model parameters are summarized in Table 7. 
Based on the value for the visual extinction Ay we interpolated the extinction for the other 
wavelengths following the findings of Mathis (1990). The disk inclination % was derived from 
the ellipses fitted to the FWHM of assumed Gaussian brightness distributions (see section 
5.3.). For the position angle 6 we used two values: One was also derived in section 5.3., and 
for the second one we fitted the disk model to the observed visibilities with 9 being the only 
free parameter (see also section 7.2.). 

It is clear that this simple approach is not able to reproduce all features of the SED 
(e.g., the weak silicate emission feature). For the moment, however, we restrict ourselves 
to this simple approach and, as mentioned above, a detailed modeling is left for further 
investigations. It is striking that rather high disk temperatures are required to fit the SED. 
In the inner 0.35 AU no dust can withstand the high temperatures >1500 K even on the disk 
surface and only the gaseous disk component can survive. It is interesting to note that the 
power law index of ql = 0.75 of the inner disk region equals exactly that what was derived in 
analytical disk models for geometrically flat, optically thick passive disks (e.g., Adams, Lada 




(6) 




(7) 
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& Shu 1987) as well as for flat, steady state, optically thick accretion disks (e.g., Pringle 
1981). Thus, from the temperature distribution alone one can not distinguish between the 
two different cases (see also Kenyon & Hartmann 1987). However, in the context of FU 
Ori the passive disk model where the disk is merely re-processing radiation from the central 
star is unlikely and only an accretion disk model matches most observations (Hartmann 
& Kenyon 1996). A nice example is the recently detected magnetic field in the innermost 
regions of the disk surrounding FU Ori which can only be explained in the picture of a 
heavily accreting disk (Donati et al. 2005). As the disk temperature is directly related to 
the mass accretion rate M in classical accretion disk models it is possible to derive the M 
required to produce the observed high temperatures. Bell et al. (1997) analyzed the effective 
disk temperature and also the midplane temperature of an accretion disk as a function of 
the mass accretion rate. Figure 1 of their paper shows that one needs M > 10~ 5 M Q /yr to 
get close to the effective temperature we find in our model. Such high values for M were also 
found by Hartmann & Kenyon (1996) and Lachaume et al. (2003). Due to the high accretion 
rate the midplane temperature of the disk in the inner «10 AU is even much higher than 
the effective disk temperature (Bell et al. 1997). 

Concerning the temperature law farther out than 3 AU the value we adopt here for ql is 
in good agreement to what can be found in the literature for isothermal flared disks (Kenyon 
& Hartmann 1987). In their model the index value diverges from the flat solution where 
q = 0.75 at a certain point in the inner disk and approaches q = 0.5 for radii much larger 
than the stellar radius. At intermediate radii the temperature index might be approximated 
by values in between these most extreme cases depending on the actual flaring of the disk. 

In summary, the temperature distribution indicates that in the innermost disk regions 
the disk luminosity is accretion dominated whereas further out it is irradiation dominated. 
This is expected for an active accretion disk since the luminosity released by accretion is a 
much steeper function of distance to the star than that by absorption of radiation from the 
central star and the hot innermost disk regions (L accretion oc r -4 , whereas L a b sorpt i n oc r~ 2 , 
roughly). 

7.2. Model visibilities 

7.2.1. MIR results 

For this simple disk model it is straightforward to compute the visibilities as a function 
of projected baseline and wavelength. Assuming that each part of the disk surface is emitting 
as a blackbody with an effective temperature T(r) the corresponding total flux and visibility 
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for an axisymmetric face-on disk is given by (see also equations (2), (3), and following in 
Malbet et al. 2005): 



with F\ being the flux, B\ the Planck function, Jo the zeroth-order Bessel function of 
the first kind, Bp the projected baseline and d the distance to the object. 

Taking into account the inclination of the disk and the position angle as given in Table 7 
we computed the corresponding visibilities for the three different baselines. The results for 
two different position angles are shown in Figure 10. In the upper panel the disk inclination 
and the position angle are derived from the ellipses fitting in section 5.3. In the lower panel 
the position angle was fitted to the observed visibilities with a least squares fit. While the 
ellipse fitting yielded a position angle of 109. Pill. 6° a least square analysis leads to a value 
of 93.4°±6.8°. From the plots in Figure 10 one can see that changing the position angle 
mostly affects the visibility of the shortest baseline while the curves for the other baselines 
change only marginally. 

Fitting the UT3-UT4 baseline proofs quite a challenge as the observed visibility is 
increasing over the 8-13//m wavelength interval while all other modeled visibilities are de- 
creasing. Especially for this baseline better fits are expected from more sophisticated disk 
models. However, given the simplicity of the current models the fits seems in reasonable 
agreement with the observations. It clearly shows that interferometric measurements can 
put additional constraints on the structure and geometry of circumstellar disks that cannot 
be derived from SED fitting alone. 

Interestingly, Malbet et al. (2005) were able to fit their NIR visibilities of FU Ori with 
a similar disk model and also with a disk and an embedded "hot spot". They concluded 
that the latter model was more likely based on statistical arguments. The major differences 
between their models and ours is that (1) they fit the SED with single power laws for the 
temperature neglecting the flux longwards of 20/im and (2) their best fits yield disk position 
angles of 47° and 8° ± 21°. The inclination angles (55° t 5 7 l and 48 l?o°» respectively), 
are in good agreement with our value. In Figure 11 we compare both of their models to our 
observations. We find that both models fit the MIDI visibilities not as good as the model 
presented here. Even more we find that our MIR interferometric measurements rule out 
their first model of a simple accretion disk. This model predicts a higher visibility for the 
UT3-UT4 baseline than for the UT2-UT3 baseline which is not observed. The model with 




(8) 
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an embedded "hot spot" in the accretion disk shows the right trend for the visibilities, but 
the values for the UT3-UT4 and UT2-UT4 baseline are higher than for our models which 
already give slightly too high values compared to the observations. The main reason for 
the higher visibilities resulting from the model applied by Malbet et al. (2005) is the single 
power law approach which provides less MIR flux compared to our models. 

7.2.2. NIR results 

We also compared our models to the measured NIR visibilities from Malbet et al. (2005). 
We restricted ourselves to those observations where the object was clearly resolved. The 
results are given in Table 8. Except for the observations done with north-west baseline 
of the Palomar Testbed Interferometer (PTI/NW) both of our models agree with the NIR 
visibilities within the error bars. 

Concerning the predicted "hot spot" in the accretion disk (Malbet et al. 2005) our 
observations do not yield new insights into its nature as they do not have the required 
spatial resolution to confirm its existence. Also the calibrated phase of our observations do 
not contain any information on which basis we could speculate about this possible second 
companion. 

8. Conclusions and Future Prospects 

We presented the first multi-baseline MIR interferometric observations of FU Orionis 
with MIDI/VLTI. The findings can be summarized as follows: 

1. FU Orionis was clearly resolved in the MIR with two VLTI baselines and marginally re- 
solved with a third baseline indicating the presence of warm dusty material surrounding 
FU Ori out to several tens of AU. 

2. The inclination and the position angle of the accretion disk can be inferred from multi- 
baseline measurements with MIDI/VLTI. Thus, the instrument provides means for 
deriving disk parameters that are otherwise only poorly constrained from models and 
SED fits. 

3. The correlated flux indicates that 95% of the 8-13/xm flux comes from within the inner 
25AU of the disk for the shortest of our baselines while 65% of this flux arise from the 
inner 13AU for our longest baseline. 
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4. The shapes and strengths of the total 8-13/xm spectrum and the (spatially resolved) 
correlated spectra indicate that most dust particles within the accretion disk are amor- 
phous and already significantly larger than typical particles in the ISM. 

5. No spectra, neither the total nor the spatially resolved correlated spectra bear signifi- 
cant traces of crystalline silicates. Given the high accretion rate and disk temperature 
of the system this is unexpected and requires further investigations. 

6. The SED and the observed MIR visibilities can be fitted reasonably well with a simple 
analytical disk model prescribing a broken power law for the effective temperature 
distribution of the protoplanetary disk. Within the innermost 3AU the temperature 
decreases with T oc r~ - 75 , farther out the temperature goes with T oc r -0 53 . The 
derived power law index for the inner regions is equal to what is derived theoretically 
for flat, steady state, optically thick accretion disks (Pringle 1981). For the outer disk 
region the derived value is in good agreement to what can be found for isothermal 
flared disks (Kenyon & Hartmann 1987). 

7. Based on points 4. and 5. and recent spectroscopic observations of protoplanetary 
disks around TTauri stars and HAeBes it seems questionable that FU Ori is a very 
young classical TTauri star as the size of amorphous dust grains in combination with 
the non-detection of crystalline dust contradicts this assumption. 

8. From the acquisition image N-band aperture photometry could be derived for the 
companion FU Ori S. It shows that this young K-type star has more relative MIR 
excess than FU Ori itself hinting at a circumstellar disk with larger flaring angle or 
smaller disk inclination. 

Putting these results into context with other high spatial resolution studies of FU Ori 
and other FUORs we see that: 

1. Most of the published NIR visibilities of FU Ori (Malbet et al. 2005) do agree with 
our newly derived disk models. 

2. Our interfero metric observations rule out one out of two disk models for FU Ori pre- 
sented by Malbet et al. (2005). Due to the lack of sufficient spatial resolution the 
presence of the proposed second very close companion can neither be confirmed nor 
disproved with the current set of MIDI data. 

3. In contrast to NIR interferometric observations of other FUORs (Millan-Gabet et al. 
2006) our data do not seem to require the introduction of an extended dusty structure 
providing additional uncorrelated flux. 
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4. For the inner 3AU we require a significantly different disk model (e.g., £ oc r~ - 9 , 
T oc r~°' 75 ) to fit our visibilities and SED in comparison to what was found for the 
outbursting star V1647 Ori (E oc r~ L5 , T oc r"°- 53 ) (Abraham et al. 2006). 

Although the present results provide important new insights into the FU Ori system 
new questions arose as well. Future investigations will include a more thorough and detailed 
modeling of the observations with numerical radiative transfer disk models. These models 
will not only need to consider the high accretion rates but also such effects as self-irradiation 
of the disk. It will be interesting to see whether these models are able to confirm the apparent 
rotation of the disk's position angle with wavelength and whether they can reproduce the 
increasing visibility curve for the UT3-UT4 baseline. Also, these disk models will have to 
reproduce the observed NIR and MIR visibilities simultaneously with higher accuracy than 
the current simple models are able to do. From the observational side re-measuring the UT3- 
UT4 baseline might be eligible and any other new baseline configuration will certainly help to 
constrain the disk geometry (especially the position angle) even better. Concerning the dust 
structure the lack of crystalline silicates needs to be analyzed in greater detail. Is it still just 
a contrast effect that we do not see any crystalline features in the spectra or are they simply 
not there? Eventually we will have to face the problem that currently all FUORs that have 
been observed with high spatial resolution techniques draw a rather inhomogeneous picture 
of the group. Apparently, some are still surrounded by nearby dense envelopes while for 
others the circumstellar disks suffice to explain the observations. However, even those disk- 
systems (i.e., FU Ori and V1647 Ori) seem to be different in terms of physical properties as 
different temperature and surface density profiles are derived. 

Without any doubt interferometric observations of circumstellar disks in the NIR and 
MIR provide an unprecedented means for deriving geometrical disk properties and put new 
constraints on the physical (and also chemical) processes taking place within protoplanetary 
disks. And as more and more interferometric data of circumstellar disks are published our 
view on the cradle of planetary systems will be once more refined. 
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Table 1. Journal of MIDI observations of FU Orionis and the calibrators used for the data 

reduction. 



Date 


Object 


Projected Baseline 


Position Angle 


Comment 


31.10.2004 


FU Ori 


44.56 m (UT2-UT3) 


46.54° 








HD37160 








Calibrator 


02.11.2004 


FU Ori 


86.25 m (UT2-UT4) 


84.24° 








HD 31421 








Calibrator 




HD 37160 








Calibrator 




HD 50778 








Calibrator 


04.11.2004 


FU Ori 


56.74 m (UT3-UT4) 


106.64° 








HD 25604 








Calibrator 




HD 20644 








Calibrator 




HD 37160 








Calibrator 




HD 50778 








Calibrator 


29.12.2004 


FU Ori 


44.80 m (UT2-UT3) 


46.61° 








HD 37160 








Calibrator 




HD 94510 








Calibrator 8 - 



a Observed twice that night. 
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Table 2. Photometric values for the FU Ori system. J to L' are taken from Reipurth & 
Aspin (2004), N-band values are derived from MIDI 8.7/xm acquisition images. The errors 
for the N-band fluxes are standard deviations based on measurements of three acquisition 
images. Two were taken on 31.10.2004 and one on 02.11.2004. 





FU Ori 


FU Ori S 


J [mag] 


6.30±0.03 


10.75±0.23 


H [mag] 


5.64±0.05 


9.92±0.21 


K' [mag] 


5.25±0.02 


9.15±0.15 


L' [mag] 


4.18±0.04 


8.09±0.16 


N [mag] 


2.75±0.19 


5.28±0.11 
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Table 3. Dust species used in the dust model. Apart from the name, the chemical 
formula, the shape and also the reference to laboratory measurements for the optical 
properties are given. Mie-theory was used to calculated the opacities for the homogeneous 
spheres whereas for the inhomogeneous spheres we used the distribution of hollow spheres 
given by Min et al. (2005) to simulate grains that are not perfectly symmetric. References: 
(1) Dorschner et al. (1995), (2) Servoin k Piriou (1973), (3) Jager et al. (1998), (4) Spitzer 

& Kleinman (1960). 



Species 


Chemical 
Formula 


Shape 


Ref 


Amorphous Olivine 


MgFeSi0 4 


Homogeneous 


(1) 


Amorphous Pyroxene 


MgFeSi 2 6 


Homogeneous 


(1) 


Crystalline Forsterite 


Mg 2 Si0 4 


Inhomogeneous 


(2) 


Crystalline Enstatite 


MgSi0 3 


Inhomogeneous 


(3) 


Amorphous Silica 


Si0 2 


Inhomogeneous 


(4) 
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Table 4. Possible dust composition for the observed 10/zm silicate feature as derived with 
our dust model. The mass fraction of three grain sizes for different dust species is given. 





0.1 /iin 


1.5 /xm 


6.0 /iin 


Total 


Amorphous Olivine 


0.17 


< 0.01 


0.18 


0.35 


Amorphous Pyroxene 


< 0.01 


0.07 


0.57 


0.64 


Crystalline Silicates a 


< 0.01 


< 0.01 


< 0.01 


< 0.01 



a Forsterite and Enstatite 
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Table 5. Calibrated visibilities for FU Ori observed with three different baselines. The 
errors are standard deviations resulting from calibrations with different calibrator stars for 

each night. 



Wavelength [/im] V [UT2-UT3 a ] V [UT2-UT4 b ] V [UT3-UT4 C ] 



1 9 


n Q«-i-n n.9 

u.yozLU.uz 


n nfH-D D9 

U.UUZLZU.UZ 


n «zi-i-n ri9 

u.o^izizu.uz 


19 8^ 




fl fil _l_fl no. 

U.U1ZLZU.UO 


fl SQ-l-fl no 
U.OOZLU.U^ 


19 7n 


pi Qzl-|-n fl^ 


f) fiQ-|-f| 


fl SzL+n D9 
u.o^iznu.uz, 


1 9 ^7 


n Qfi-L-f| D9 
u.yuzLU.uz 


n R9-|-n no 

U.U-^ZLU.U.^ 


n szL+n ri9 

u.o^zlu.u^ 


1 9 zL/L 




n n 9+n n9 

U.UZZLU.UZ 


fl 01 4_n no 
u.oizlu.uz 


1 9 ^1 


n QzL-i-n n9 


u.uuzcu.ui 


fl Ol _L.fl fi9 

u.oiznu.uz 


19 17 


u. c/z,zizu .u 1 


n «i _i_n no 
u .u izcu. 


fl 7q-L.f1 fl9 

U. 1 C/ZLZU.UZ, 


19 03 


D QzL+n 09 


n fil -4-fl 09 
U.uizcu.uz 


fl 01 _L.fl fi9 

U.OIZLZU.UZ 


1 1 


n QQ_i_n no 
u.yozizu.uo 


n fi94-f| 09 
u.u^zlzu.uz 


fl 7Q-L-fl fl9 

u. i c/ZLzu.uz 


1 1 7^ 
11.1 o 


n QzL+n 


n n O-|-n ni 


fl 7S-L-fi fi9 

U. ( OZLZU.U^ 


1 1 fin 


u. yixu.uo 


n n Q-i-n ni 

u.uoznu.ui 


fl 774-fi fl9 

U. 1 1 ZLZU.UZ, 


1 1 zL^ 


n 01 _|_n f|9 
u.c/1zlzu.uz 


U.UOZLZU.U1 


n 7Q-L-fi fi9 
u. / yztzu.u^ 




n oa _i_n riQ 


n «9-un D1 
u.uuznu.ui 


fl 7S-L.fi f)9 
u. t oztzu.uz, 


1 1 1 A 


n QO-|-n 
u.y^zizu.uo 


u.uzznu.ui 


fl 7fi-L.f1 fi9 

U. 1 UZLZU.U^ 


1 n q& 


n QQ-L.fi n9 
u.yozizu.uz 


n «o_i_n ni 
U.uuznu.ui 


fl 774-fi fi9 
u. 1 1 zLzu.u^ 


1 n &9 


n Q9_un n9 
u.yznzu.uz 


fl fiO-Lfl fi1 
U.UZZlU.UI 


fl 7fi-L.fl ni 
U. ( UZLZU.Ul 


10.66 


/ 1 no 1 n no 

0.92±0.02 


n /-o i n n 1 

0.62±0.01 


n *rr 1 n ni 

0.75±0.01 


10.49 


0.91±0.01 


0.62±0.01 


0.75±0.02 


10.32 


0.92±0.02 


0.63±0.01 


0.75±0.02 


10.15 


0.92±0.01 


0.63±0.01 


0.75±0.02 


9.97 


0.94±0.02 


0.65±0.01 


0.73±0.01 


9.80 


0.97±0.02 


0.64±0.01 


0.73±0.01 


9.62 


0.97±0.03 


0.64±0.01 


0.74±0.01 


9.43 


0.93±0.02 


0.61±0.01 


0.76±0.02 


9.25 


0.93±0.01 


0.67±0.01 


0.76±0.02 


9.06 


0.93±0.01 


0.69±0.01 


0.74±0.02 


8.86 


0.94±0.01 


0.69±0.01 


0.74±0.02 


8.67 


0.95±0.01 


0.70±0.01 


0.73±0.02 


8.47 


0.95±0.01 


0.70±0.01 


0.74±0.02 


8.27 


0.94±0.01 


0.73±0.01 


0.73±0.01 
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Table 5 — Continued 



Wavelength [/im] V [UT2-UT3 a ] V [UT2-UT4 b ] V [UT3-UT4 C ] 



a Projected Baseline: 44.56 m, Position Angle: 46.54° 
b Projected Baseline: 86.25 m, Position Angle: 84.24° 
c Projected Baseline: 56.74 m, Position Angle: 106.64° 
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Table 6. Upper half: Derived FWHM (in AU) of Gaussian brightness distributions for the 
three baselines and three wavelengths using equation 2. The distance to FU Ori is assumed 
to be 450 pc. Lower half: Parameters of best ellipses fitted to the data for the three 

considered wavelengths. 



9.0 jum 11.0 jum 12.5 /im 



UT2-UT3 (44.56 m) 
UT2-UT4 (86.25 m) 
UT3-UT4 (56.74 m) 


o 47+0.40 

-0.46 

O OC+0.24 

°- zd -0.25 
4 1 O+0.28 
^• iy -0.30 


q 97+0.67 
' -0.82 

4.34^;° 8 7 
5.0018$ 


o 40+O.8I 
°-^°-1.02 

4-9318;;? 

4.83i8:g 


Semimajor axis a [AU] 


2.28 


2.51 


2.54 


Semiminor axis b [AU] 


1.21 


1.54 


1.42 


Ellipse area [AU 2 ] 


8.67 


12.14 


11.33 


Position angle # a 


122.2° 


111.1° 


93.9° 


Inclination i° 


58.0° 


52.2° 


56.0° 



a Position of a measured from north eastwards 
b Assuming an underlying circular disk 
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Table 7. Parameters of the simple disk model used to fit the SED. The upper part of the 
table gives the parameters for the inner 3 AU. Farther out in the disk the temperature and 
surface density follow distributions defined by the parameters in the lower part of the 

table. 



Parameter 


Variable 


Value 


Inner disk radius 


Rin 


5.5 R Q 


Outer disk radius 


Rout 


100 AU 


Disk inclination a 


i 


55.4° 


Disk position angle b 


9 


109.1° a / 93.4° c 


Extinction 


Ay 


2.6 


Temperature at 1 AU 


T\AU,\ 


670 K 


Power law index for temperature 


ql 


-0.75 


Power law index for surface density 


pi 


-0.9 


Disk mass 


M D>1 


0.03 M 


Temperature at 1 AU 


T\AU,2 


550 K 


Power law index for temperature 


q2 


-0.53 


Power law index for surface density 


p2 


-1.4 


Disk mass 


M Dj2 


0.01 M 



a As derived in section 5.3. 

b Only required for visibility computation 

c Resulting from x 2 ~fit 
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Table 8. Comparison of observed NIR square visibilities from Malbet et al. (2005) to our 
disk model with two different position angles. The observations were carried out with 
different baselines at the Palomar Testbed Interferometer (PTI) and the Very Large 

Telescope Interferometer (VLTI). 



Interferometer 


Filter 


Baseline [m] 


Pos. Angle [deg] 


V 2 (observed) 


V 2 (model) a 


PTI/NS 


H 


103.7 


61.8 


0.83±0.04 


0.80 / 0.82 


PTI/NW 


H 


85.6 


14.1 


0.79±0.05 


0.92 / 0.94 


PTI/NS 


K 


102.7 


63.6 


0.72±0.08 


0.78 / 0.80 


PTI/NW 


K 


84.2 


14.0 


0.79±0.05 


0.91 / 0.93 


PTI/SW 


K 


82.5 


-55.6 


0.82±0.08 


0.85 / 0.87 


VLTI/UTl-UT3 b 


K 


89.6 


54.2 


0.87±0.05 


0.83 / 0.85 



a For two different position angles: 109.1° / 93.4° 
b VINCI data 



31 




Fig. 1. — FU Ori and FU Ori S as seen in the MIDI acquisition images at 8.7//m. The image 
was taken with UT3 the 31st of October 2004 using the AO-system MACAO (north is up, 
east to the left). 
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Wavelength [urn] 



Fig. 2. — MIR spectrum of FU Ori from 5.3-14/iin. Data from the MIDI observations (dots 
with solid line) range from 8-13/im. They agree within the error bars with Spitzer IRS data 
(squares). The errors for Spitzer are smaller than the symbols and are not shown. The larger 
MIDI errors arise from averaging six independently calibrated measurements (3 nights and 
2 telescopes each night). The observed silicate feature is rather flat and broad. 
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Fig. 3. — MIR spectrum of FU Ori and the resulting dust model fit. The Spitzer spectrum 
(black solid line) can be fitted quite well with the dust composition given in Table 4. The 
dotted line shows the contribution of the underlying continuum, the dash-dotted line shows 
the silicate emission and the grey dashed line shows the sum of both, i.e. the total flux, 
which is hardly distinguishable from the observed spectrum. 
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o-e UT2-UT3: 44.56 m, 46.54 deg 
q- -a UT2-UT4: 86.25 m, 84.24 deg 
UT3-UT4: 56.74 m, 106.64 deg 



8.5 9 9.5 10 10.5 11 11.5 12 12.5 13 
Wavelength [|xm] 



Fig. 4. — Calibrated Visibilities of FU Ori for three different baselines. 
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Fig. 5. — Sizes and orientation of the emitting regions on the sky for different wavelengths 
(9.0yum (upper left), 11.0/xm (upper right) and 12.5//m (lower left)). The sizes correspond to 
FWHM of assumed Gaussian brightness distributions and are given in Table 6. The resulting 
best fit ellipses for each wavelength are overplotted and combined in a fourth plot (lower 
right). The center of the plots corresponds to the position of FU Ori and north is up and 
east to the left. 
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Fig. 6- 



Observed correlated flux for the three different baselines and position angles. 
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Fig. 7. — Upper Panel: Continuum subtracted and normalized correlated spectra for the 
three different baselines and the normalized total spectrum. For clarity the errorbars have not 
been overplotted. They equal those of Figure 6. Lower Panel: Variations of the normalized 
correlated spectra with respect to the normalized total Spitzer spectrum in units of error. 
It shows that most variations are clearly below the 3a threshold, indicating that the shape 
of the correlated spectra equals that of the total spectrum. 
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Fig. 8. — Spitzer MIR spectrum from 14-34/xm. The spectrum was smoothed by a factor 
of three. The typical errors are smaller than the size of the symbols. Similar to the shorter 
wavelength spectrum in Figure 2 this spectrum appears rather smooth without any striking 
features. The arrows indicate where the order of the spectrograph is changing leading to 
variations in the detected flux (approximately at 19.6, 22.1, 25.0, and 26.7//m). Although the 
weak peak around 27.7-27.8/im can theoretically be attributed to the presence of forsterite 
particles their detection seems questionable as the normally stronger emission bands at 23.9 
and 33.8/^m is clearly not seen. 




-39- 



le-10 

le-11 

le-12 

S le-13 
I 

Ph le-14 
le-15 
le-16 



iV ^' 1 10 100 1000 

Wavelength [|Xm] 

Fig. 9. — SED fit of a geometrically thin, optically thick accretion disk model. The disk 
parameters are shown in Table 7. The contribution of the central star (in this case we assume 
T e ff = 3700 K, R s tar = 2.5R Q ) can be neglected. Although simple, the disk model agrees 
reasonably well with the observed fluxes for most wavelength regimes. The data points were 
taken from Clarke et al. (2005) (U,B,V,R), Przygodda (2004) (8.9/mi, 11.9/mi), Reipurth & 
Aspin (2004) (J, H, K, L), Weintraub et al. (1991) (450/rni, 850/im, 1.3mm), Spitzer IRS 
data archive (MIR Spectrum) and IRAS data archive (12 /mi, 25 /mi, 60 /mi, 100 /mi). 
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Assumed Disk Position Angle: 109.1 deg 
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Assumed Disk Position Angle: 93.4 deg 
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Fig. 10. — MIR visibilities derived from the simple disk model in comparison to the obser- 
vations. The baselines and projection angles correspond to those given in Table 1. The disk 
position angle was assumed to be 109.1° in the upper panel and 93.4° in the lower panel. 
Changing the position angle has the biggest impact on the visibilities of the shortest baseline. 
Apparently, both models do show the right trend but they can not account for increasing 
visibilities. 
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Fig. 11. — MIR visibilities derived from models for FU Ori presented in Malbet et al. (2005) 
in comparison to the observations. Upper Panel: Model consisting of a simple accretion 
disk. This model can be ruled out by our observations as it predicts higher visibilities for 
the UT3-UT4 baseline than for the UT2-UT3 baseline which is not observed. Lower Panel: 
Model consisting of an accretion disk and an embedded "hot spot". This model shows the 
right trend but does not fit the observations as good as the models presented in this paper. 



